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Abstract 

The stochastic vacuum model description of a heavy meson is discussed 
in the context of a gauge-invariant approach where Wilson loop expectation 
values appear naturally in the 0(v 2 ) spin-orbit Hamiltonian. These expec- 
tation values have been derived elsewhere, however by a procedure whose 
legitimacy is now placed in question. Here they are derived by standard func- 
tional methods with a result that is identical to the previous one. In addition, 
a full spin- independent Hamiltonian reduction to 0(v 2 ) is carried out. 



1 introduction 



It was the work of Leutwyler and Voloshin some years ago in the context of the 
sum rule formalism that first suggested the fundamental nonlocality of nonpertur- 
bative interactions between hadronic quarks 0. Leading effects were later shown 
proportional to a gluon condensate thereby excluding the possibility of a purely 
local description. In a separate line of development Wilson's lattice work || led to 
the well-known area law as a qualitative formulation of color confinement. There, 
gauge invariance of the hadronic state is the guiding principle. There too what is 
fundamental is the presence of a nonlocal operator - the Wilson loop. 

In the stochastic vacuum model (svm) of Dosch and Simonov M we find these 
salient features mutually complementary, where area law asymptotics follow from a 
non-zero gluon condensate. Active gluon degrees via the field's correlation length 
measured against QQ correlations also come into the picture; this in such a way 
that the model also accounts for intermediate as well as short-range perturbative 
behavior. 

The aim of the present article is to describe the reduction of the svm to an 
effective 0(v 2 ) spin-orbit interaction Hamiltonian for the heavy QQ state. In fact, 
a general reduction in terms of Wilson loop expectation values (ev) has already 
been made ||, thus leaving as the main focus of this study only the ev derivations 
themselves. The ev too have been derived [H 0, H P> 0) however by a procedure 
whose legitimacy is here placed in question: From the functional variation of the 
Wilson-loop 

5i\nW = J 5a^( Zl )((F^(zi))) (1) 
is presumed the functional derivative relation 

resulting in a spin-orbit interaction that e.g. satisfies the Gromes relation. 
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We argue that procedure (1) — > (2) lacks legitimacy by virtue of insufficient 
degrees of freedom characterizing the variational area element 5a llll . In section 2 the 
familiar svm spin-orbit Hamiltonian is derived along the usual lines from assumption 
(1) — > (2). We then demonstrate the general invalidity of the assumption and show 
that when omitted it leaves the Hamiltonian insufficiently specified. Additional 
constraints are derived in section 3 where the model statement itself is expanded to 
0(m~ 2 ) and velocity coefficients matched term by term. This procedure leads again 
to the familiar svm spin-orbit Hamiltonian; I.e, our spin-orbit result agrees with the 
one derived from (2), although, as we point out, generally applied (2) may lead to 
incorrect results [O]. Finally, we lay out details of a full reduction to the 0(m~ 2 ) 
spin-independent Hamiltonian. 



2 functional variation and the svm 

From a Foldy-Wouthysen reduction of the gauge invariant QQ 4-point function the 
0(m~ 2 ) interaction Hamiltonian of || is given in terms of the Wilson loop 

V = V so + V si 

Vso = E(-i) i ^-^ feS ^r((F^(^)}) 

1 ' Lx • s x - -^L 2 • [V (r) + 2V 1 (r)]'/r 



2ml 2m 2 
1 



(Li • s 2 - L 2 • s 1 )Vi(r)/r (3) 



J dtV si = ihiW = iln-(trPexp(igfdt(Ao-z i A i ))) 

where Darwin and hyperfine terms, not relevant to the present discussion, are omit- 
ted. Key to this approach therefore is the evaluation of the six independent expec- 
tation values ((F^ v )). The defining Euclidean svm statement in the present context 
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IS 



-~ / da pv {u) I da Xp (v)\(S pX 5 up - Sp P S uX )D(w 2 ) 



(4) 



d d 
-^-{w x 5 up - w p 5 vX ) + ^(wpd^ - w x 8 w ) 



Z^i(w 2 )|, w = u — v 



where the integrals are evaluated over instantaneous straight-line surfaces: (u{ = 
szu + (1 — s)z 2 i,u 4: — t), < s < 1, with da pv {u) = dtds(du fl / dt) {du u / ds) . f3 is the 
gluon condensate, and D and D\ are gluon correlation functions that fall off rapidly. 
From the quark world-line variation, of the Wilson loop 



8i\nW 



fcv(*iK<*W(*i)>> 
-if3 / Sa^izt) / da Xp (v)\(5 pX 5 up - 5 lip 5 uX )D(w 2 ) 



(5) 



dz 



d „ d 

{w x 6 up - w p 5 uX ) + - — [w p S pX - w x S pp ) 



1(« 



dz u 



W = Zi — V 



is proposed in references || [7|, |8| the following the functional relation 
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-i\nW 



(6) 



where 5a ^ = [dz^bzy — dz u Sz ll )/2 is called the variational area element. From this 
follows 



((Faizi))) = fa I dr{^d\ l -D{T 2 + \ 2 )+ l -D l {r 2 + r 2 ) 



(7) 



{{F u {zx))) = (3(zuri - z u ri) I dr [ d\- 

Jo r 



D(t 2 + A 2 



A 



+(3(z 2l r t - z 2l n) dr{ d\-D(r 2 + A 2 ) + -D^r 2 + r 2 ) 



o r L 



yielding spin-orbit potentials 

Vftr) = pfdr 



rfA J D(r 2 + A 2 ) + -D 1 (r 2 + r 2 ^ 
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V((r) = -(3 J dr jT d\(^l-^jD(T 2 + A 2 ) 
Vi(r) = pfdr 



(8) 



dA ^>(T 2 + A 2 ) + ~L> 1 (^ 2 + r 2 ) 
or 2 



which together e.g. satisfy the relation of Gromes [12 



[Vv + Vx-Vz]' = 0. 



(9) 



We show here that variational derivative (||) does indeed not follow from variation 
(|5|). One should notice that in passing from (|5|) to (|6]), (1) to (2), is the assumption 
that there are in the variational area element at least six quark and six antiquark 
coordinate degrees of freedom - an assumption clearly in error. 

The Wilson loop is a functional of quark and antiquark world lines. With this 
in mind we rewrite the rhs of fl5|) 



dtdzxvZx^F^zi))) = J dtSz^z^Q^ 



(10) 



where 







r.2\ 



-i/3 j dcr X p(v)\(5p\5 vp - 5 w 5 v \)D(w 
1 [ d „ d 



DAw 2 ) 



From this we should like to extract the six element {{(F^))} for determination of 
the Hamiltonian (|3|). The set appears in (|^) as the linearly transformed 



{SijSp^p ( (F uj ( Zl ))), 6 ijfl vZin ( (Fvj{zx )))} 



(11) 



- a set of vectors spanning the space of field tensor elements. The first subset 
accounts for spin independent V si , by Stokes theorem, and the second for spin de- 
pendent V so directly. It should be clear that the first subset appears in fl5|). Hence in- 
dependence of V so with respect to variation (151) follows from the linear independence 
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of vectors in the complete set (|11|). The independence is immediately established 
from the nonsingularity of the set's coefficient matrix 



det 



V 



1 








-z 2 


-z 3 








1 





Zl 





-z 3 








1 





Zl 


Zl 


-z 2 


k 





2 








-Z3 





Zl 





2 








-z 3 


z 2 








2 



~2\2 



(2 - 



where the vectors have been arranged by row with field tensor elements in ascending 
order, left to right. 

For concreteness, then, variation (|5]) 



dt5zuZi j ((F4 j (zi))) 



dtSzuZij&4j (12) 
J dt5z lj (Q 4:j + zuQij) 



leads to 



«*Q,-(*l)» 

((Fijizi))) 



®oj + {zfrj - zijZiiri)fi{r) + {zuz^rj - z 2 jZuri)f 2 {r) 



<S>ij + (hzij - fjZii)f[(r) + {fiz 2 j - rjZ2i)f 2 {r) 
yielding when installed into (Q) the spin-orbit potentials 

V((r) = f[{r)-(3 j dr j^d\(l-^D{T 2 + \ 2 ) 



(13) 



(14) 



V'(r) 



-f 2 {r)+0 dr 



dX ^D(T 2 + \ 2 ) + r -D 1 (T 2 + r 2 ) 
r 2 



for arbitrary functions (f[,f 2 )- I.e., solutions (|H]|) and (TJ) satisfy variation (|5|), 
yet leaves the spin-orbit Hamiltonian V so entirely unspecified. This claim may be 
verified by direct substitution of (|i~3"D into (|)and (p~3|) into @. 
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3 specification of the svm Wilson loop expecta- 
tion values 



We now revisit the spin-orbit Hamiltonian derivation beginning again from (0) . As in 
ref ||1 1|| , additional constraints are found from the functional expansion of the model's 
defining statement in orders of heavy quark velocity. Accordingly, definition @ is 
expanded 

i In W ~ J dtdt'{v°(r, r') + (z u + z 2i )V°(r, r') (15) 
+ + z 2i z' 2j )V?(r,r') + { Zli z' 2j + z 2i z' Xj )Vri{r, r')}. 

which on the 0(v 2 ) approximation 

r\ « n - fiT + fir 2 /2 , T = t-t' (16) 

yields 

duty ~ | df|y (r) + (z? + z 2 2 )V a (r) + z, ■ z 2 V b (r) (17) 
+ [( Zl • rf + (z 2 • f) 2 ]V;(r) + (zi • f)(z 2 • r)V d (r)}. 

for functions y a , V^ 3 , \^ 7 , Vo, H, K> ^ to be determined. Equation ( [15] ) leads to 
( see appendix ) 



{{F^zx))) = (3(z lj r i -z li r J )JdrJ^dX-i [ l--jD(r 2 + X 2 ) (18) 
+P(z ij r i - z 2l r 3 ) J dr ij^dX^D(T 2 + A 2 ) + l -D x {r 2 + r 2 )| , 



and for the spin-independent potentials of (|17D we find upon application of ([R]) (see 
appendix) 



V (r) = 13 dr dX 



o 



(r-X)D(r 2 + X 2 ) + ^D 1 (r 2 + X 2 ) 



(19) 
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V a (r) 



-P / dr 



Vh(r) 



dX 

fx 



-p dr dX 
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Ar 2 


L \b 


4 




0^2 


X 2 


A 3 


Ar 2 


r 2 


4r 


4r 2 


4r 


- / r 


A 3 


Ar 2 






6r 2 




r y 


X 2 


A 3 


Ar 2 


-V 


2r 


4r 2 


2r 2 


2r / 



(20) 



D(r 2 + A 2 ) 



VJr) 



VJr) 



~P 

-P I dr 



dr / dX 



r Xr_ _ Xr_ 
6 2r 3r 2 



dA 



r 



A 2 2A 3 \ 



-- — + — ) Dir + a ) - + r 



(21) 

A 2 )l (22) 

6 r ' 3r 2 / ~ v ' ' " y 4~ 1V ' ' y J ^ 
A couple observations: The ev of fll8|) follow from strict independence of operators 
(r, f) in the order expansion (15). The independence is relaxed in the 0(t> 2 ) evalu- 
ation (0) that leads to spin- independent potentials (|19D - (p3| ). This explains the 
discrepancy between ev flI8|) and the one obtained in an earlier effort by the present 
author |TB[ where instead the ev are derived improperly from expansion QT7p. We 
note that as (|T8|) is identical with fl7|) it yields upon insertion into (§) the spin-orbit 
potentials (§) that satisfy relation @. However, only incidentally fl3fl . We also note 
that in the long-range regime the spin-independent interaction, ( |I~9"D - (0), agrees 
with minimum area or flux-tube asymtotics || |TT|] 

1 



1 L 2 
•-a — 

6 r 



1 1 
+ — 



2 

mf 



17117712 



(24) 



a = /3y dr dXD(r 2 + X 2 ) . 

This too is a correction to the aforementioned earlier effort where in the correspond- 
ing equation ([17]) the relation 

7 2 = -7 2 (9*>) 

Euclidean Minkowskian \ u ! 

is overlooked. In ref.|§ a spin independent Hamiltonian reduction of the svm to 
0{v 2 ) is carried out. The apparent disagreement with potentials ([19]) - (|23|) is due 



S 



to the use of a partial integration scheme slightly different from the one applied 
here. The two results are equivalent. 



4 appendix 

From the functional identities 



S(t - f) (26) 



= ^(|/(0) (27) 



the Taylor function and functional expansions are carried out to 0(v 2 ) 

/(Z!,z 2 ) = /(0) + ij / df ^/(z;,z 2 )^) df QL/( Zl ,z^ 



2 L L J J \5z\5z'i 

+z\zi J J dt'dt" (^|^jM^2)) (28) 
F[ Zl ,z 2 ] = F[0]+ JdtzlfJ-Flz'^z^ +Jdtzi(J-F[z 1 ,z' 2 ^ 

+i//^if(^K,*j) o 



27 7 ^ 1^2^ /o 
+ //^^(^F[*i^) o+ ^. (29) 

where primes indicate time dependence, e.g., £' = £(£') , and subscript " 0" means 
evaluation at zq = z 2 = 0. The spatial field tensor Wilson loop expectation value to 
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first order is then 



5 



((Fifa))) = «i^» + zik / dt' — + z 2k dt' — — ((Fij))' 



6z lk 

((^■))o-((^A))o 

5 , , <) 



+i 2 , / eft' ( r -|-<(fl^»i - t|-((^}}' 1 



(30) 



To find the above rhs we expand both sides of the svm, (||), to 0(t> 2 ) 
i\nW ~ zlnWo- /dt(iu((A(^i)}}o + i2i((A(^))}o) 



5 



5 



5 



5 



Sz 



Li 





cftcft'|\/ (r, r') + (i H + i 2i )Vf (r, r') 
+(z u z' Xj + z 2i z' 2j )V?(r,r') + (z u z 2j + z 2i z l:j )V^(r,r') 



bz- 



1 1 



where we find 



f 1 , d 

V* = —i/3 J dsds's^iTjr'j — Ujrjr^r-^-^Di 

= ^ i 2j Q dsds ss U5ijr k r k - r-r^D + D t ] 



V; 



i.i 







x \uJkUirkr\5ij + u^jr^ ~ UiUJ k rjr k - u j Ukr k r' i \—D 1 



-d £dsds'(l - s)s'\{5 l3 r k r' k - r[r 3 ){D + D x ) 



(31) 

(32) 
(33) 

(34) 







x \ukUJirkr\5ij + LdiUJjr k r k - UiU k rjr' k - u j Ukr k r' i \—D 1 



and where 



0Z\j 0Z\j 



—ilnl(trPexjp[ig [ dt'{A Q + i^)]) 

2i i O J 
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= -((M*i))) (35) 
has been used. From this follows (equating velocity coefficients) 

«M*i)))o = - J dt'V?(r,r') 

dt \jh {{AjiZl))Y ) = - 2 1 dt ' v ^ r ') ■ ( 36 ) 

This result installed into (|30|) yields the ev (|TS[) upon performing the indicated 



differentiations. 

For the spin- independent Hamiltonian we make the further 0(t> 2 ) reduction on 
the rhs of (gD 

r\ ~ Ti — fjT + fjT 2 /2 

zV) - (i + o4 + °«^^ 2 ) fl ( r! + AV » (37 > 



where 



O i:j = Ar(r i Uj + r^) + t 2 ^ - Ar 2 (r i i) j + r i i) i )/2 

Uj = S2^ + (1 - s)z 2 i 

fi = s'zii + (1 - s')z 2 i 

A = s - s' , t = t - t' (38) 



which leads from (RTI) to the spin-independent potentials ([19]) - (|23|) upon performing 



the necessary summations and partial integrations. 
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